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1. Summary 

To identify near complete sets of genes required for the cell cycle and cell 
shape, we have visually screened a genome-wide gene deletion library of 
4843 fission yeast deletion mutants (95.7% of total protein encoding genes) 
for their effects on these processes. A total of 513 genes have been identified 
as being required for cell cycle progression, 276 of which have not been pre- 
viously described as cell cycle genes. Deletions of a further 333 genes lead to 
specific alterations in cell shape and another 524 genes result in generally mis- 
shapen cells. Here, we provide the first eukaryotic resource of gene deletions, 
which describes a near genome-wide set of genes required for the cell cycle 
and cell shape. 
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study. 

Electronic supplementary material is available 
at http://dx.doi.org/10.1098/rsob.130053. 



2. Introduction 

Understanding how cells reproduce and how they generate their shape are two 
major goals in eukaryotic cell biology. These two processes are related because, 
during the cell cycle, cells duplicate cellular components and reproduce their 
cell structure in space to generate two daughter cells. A key function of the 
cell cycle is to ensure accurate replication and segregation of the genome, 
because errors in genetic transmission can cause mutations and chromosomal 
rearrangements that may lead to cell death or disease. Failure to accurately 
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reproduce and maintain cell shape can disrupt tissue archi- 
tecture or influence cell motility and may also lead to cell 
death or disease. Given the importance of these two processes 
for cell biology, we have generated a genome-wide resource, 
cataloguing the genes that when deleted disrupt the cell cycle 
or cell shape in the fission yeast Schizosaccharomyces pombe. 
This is the first such resource that qualitatively describes a 
near complete set of genes required for these processes in a 
eukaryotic organism. 

Fission yeast is very amenable for investigating the cell 
cycle and cell shape [1-3] and has been used extensively 
for cell cycle studies for many years [4]. It is a rod-shaped, 
unicellular eukaryote that grows by apical extension and 
divides by medial fission and septation. This regular cell 
shape has made fission yeast a very useful organism to ident- 
ify genes involved in the cell cycle and the generation and 
maintenance of cell shape. Mutants can easily be identified 
by visually screening for cells that divide at a longer or 
shorter length compared with wild-type (WT) (cell cycle 
defect), or that do not have a rod-shape (cell shape defect). 
A long cell phenotype is generated if cells are blocked or 
delayed in cell cycle progression because they continue to 
grow but fail to divide and thus become elongated. However, 
not all genes required for the cell cycle show a long cell 
phenotype when deleted. For example, genes encoding 
checkpoint proteins that are not required during a normal 
cell cycle may have a WT deletion phenotype and genes 
required for mitosis often arrest with a more irregular, non- 
elongated cell shape [5]. The long cell phenotype is easily 
identified by visual screening and is definitive for cells 
blocked or delayed in cell cycle progression through the Gl, 
S, G2 and cytokinesis phases of the cell cycle [6]. We have 
therefore focused on identifying all genes with this pheno- 
type when deleted, to determine the majority of genes 
required for progression through interphase or cytokinesis 
in fission yeast. 

Genes important for the generation of cell shape are also 
easily identified, because mutants deleted for these genes lose 
the normal rod-shape. These cells exhibit a range of pheno- 
types, including rounded or stubby, curved, branched, 
skittle-shaped or more generally misshapen [2,3]. The pene- 
trance of these mutant phenotypes can be quite variable; in 
some cases, the majority of cells have the same altered 
cell shape, whereas in others the phenotype is of lower 
penetrance with fewer of the cells exhibiting the phenotype. 

Fission yeast currently has 5059 annotated protein coding 
genes, and a genome-wide deletion collection has been con- 
structed with 4836 genes deleted [7]. In this study, we have 
systematically visually screened the deletion collection to 
identify genes required during interphase of the cell cycle, 
cytokinesis and for cell shape. This resource complements 
and extends earlier studies using gene deletions in budding 
yeast [8,9] and RNAi-based gene knockdowns in metazoa 
[10-15]. Although the budding yeast gene deletion collection 
has been extensively investigated, it has not been subjected 
to a systematic screen for cell cycle genes such as that carried 
out here, while metazoan RNAi cell cycle studies have 
shown only limited reproducibility [12]. Given that 3397 
(67.14%) of the fission yeast protein coding genes have 
identifiable orthologues in metazoa (http://www.pombase. 
org), the genome-wide resource provided here will help to 
identify cell cycle and cell shape genes in other eukaryotes, 
including humans. 



3. Results 

3.1. Screening of the haploid deletion mutants 

We have microscopically examined and described the deletion 
phenotypes of 4843 haploid gene deletion mutants of both 
essential and non-essential genes, after sporulating diploid het- 
erozygous deletion mutants and germinating haploid spores 
on rich medium plates (see §5.1 for details of the screen and 
electronic supplementary material 1, tables SI (column G)- 
S3). Mutants were classified to one of 11 cell shape phenotypes 
together with three additional categories, namely: WT; arrested 
as normal spores (spores); or arrested as normal germinated 
spores (germination) (figures 1 and 2a and table 1; electronic 
supplementary material 1, table SI, column H), making 14 phe- 
notype categories in total (see §5.1). These phenotypes were 
mapped to terms in the Fission Yeast Phenotype Ontology 
(http:/ / www.obofoundry.org/ cgi-bin/ detail. cgi?id=fypo) 
(see the electronic supplementary material 1, table S4), and the 
gene list for each category can be found in the electronic sup- 
plementary material 1, table S5a-n. All 4843 deletions were 
assigned to only one of these 14 categories for analysis. The 
variation in the phenotype for mutants in a particular category 
can be found in the phenotype description in the electronic 
supplementary material 1, table SI, column G. For example, 
some long mutants may be slightly curved, and T-shaped 
cells were observed in a subset of curved cells. 

Gene Ontology (GO) term enrichment of biological pro- 
cesses and cellular components for these 14 phenotype 
categories is shown in tables 2-4 (see the electronic sup- 
plementary material 1, table S6a-n and table S7a-n for the 
complete results). There are 643/4843 genes with no GO pro- 
cess annotation. Of these 643 'unknowns', 574 (89.2%) have 
a WT deletion phenotype. This means that most genes show- 
ing one or more of the 13 other deletion phenotypes are 
assigned a biological process either by inference from other 
organisms or because they have been partially characterized 
in fission yeast. However, their cellular shape is often not 
part of that characterization. 

To identify genes required for the cell cycle and cell 
shape, we focused on the 11 cell shape categories. These phe- 
notypes and their relevance for the cell cycle and cell shape 
are described in the following sections. To demonstrate the 
use of this resource, we describe a more detailed cytological 
analysis of a subset of mutants altered in cell shape. In 
addition, we screened the library of haploid viable deletion 
mutants for hydroxyurea (HU) sensitivity and identified 
new genes implicated in the DNA checkpoint preventing 
entry into mitosis when DNA replication is incomplete or 
DNA is damaged. 

3.2. Gene deletion mutants with an elongated 
cell phenotype 

A long cell phenotype identifies cells blocked in progression 
through interphase of the cell cycle or cytokinesis (see §2). 
Gene deletion mutants with this phenotype were assigned 
to three categories, long high penetrance (long HP; 346/ 
4843), long low penetrance (long LP; 136/4843) and long 
branched (long Br; 31/4843) (figure 1 and table 1 and elec- 
tronic supplementary material 1, tables S4 and S5g-i). 
These three categories totalled 513 genes required for cell 
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Figure 1. Cell shape categories. Examples of the 11 major cell shape phenotype categories described in this study. For each deletion, mutant cells are shown as 
observed during the screen. Definitions of each phenotype can be found in the electronic supplementary material 1, table S4. Scale bar (shown in WT) = 10 [xm. 



cycle progression, and we conclude that 10 per cent of fission 
yeast genes are required directly or indirectly for progression 
through interphase or cytokinesis (513/4843). Four hundred 
and sixty-seven genes had a strong elongated deletion pheno- 
type and a further 46 genes were included that were only 
weakly elongated when deleted (see §5.1). Of the 513 genes 
in the long categories, 66.5 per cent (341 genes) were essential 
for viability compared with 26.2 per cent for all genes, and 
85.57 per cent (439 genes) were conserved in human. This 
indicates that many of the genes identified in this study are 
likely to be important for understanding the cell cycle in 
other more complex eukaryotes. 

All three long categories were enriched for nuclear localiz- 
ation (table 3). The long Br set was enriched for cytokinesis 
and transcription and included genes encoding subunits of 
the RNA polymerase II holoenzyme, mediator and SAGA 
complexes (tables 2 and 4). As expected the long HP and long 
LP sets were both overrepresented for genes involved in 
DNA metabolism and the regulation of mitotic cell cycle 
(table 2). The long HP set was also enriched for processes and 
complexes involved in mRNA metabolism (particularly spli- 
cing), RNA biogenesis, transcription and DNA replication 
initiation, and 6/6 genes encoding subunits of the MCM com- 
plex (see table 4 and legend for details). By contrast, the long LP 
set was enriched for chromosome segregation and for genes 
encoding subunits of kinetochore complexes, including Mis6- 
SIM4 (8/14 subunits), Ndc-Mis-Spc (6/10 subunits), condensin 
(4/5 subunits) and AFC (8/13 subunits) (tables 2 and 4). These 
differences in enrichment indicate that long HP is a good 



classifier for genes involved in progress through interphase, 
whereas long LP is more specific for genes associated with pro- 
gress through mitosis (table 2 footnotes g,h). Mutants that 
block in mitosis usually display an irregular less defined cell 
shape but are not elongated, whereas an elongated phenotype 
is characteristic of an interphase block, so one possibility is that 
the long LP gene set is enriched for a subset of mitotic genes that 
are also required for interphase progression, with some cells 
arresting in mitosis and other cells in interphase. Together, 
these three long phenotype categories (long HP, long LP and 
long Br) define gene sets important for progression through 
interphase or cytokinesis. The long LP gene set is in addition 
enriched for a set of mitotic genes, which may also be required 
during interphase. 

3.3. Cell cycle mutants without a long phenotype 

Not all genes previously known to be required during inter- 
phase have an elongated phenotype when deleted. We found 
that DNA replication genes cdc23, ssbl, poll, cdcl8, cdtl, rad4 
and four genes encoding core subunits of the RFC {rfc2, rfc3, 
rfc4 and rfc5) were all annotated to the misshapen essential 
set {miss E; table 1 and electronic supplementary material 1, 
table S5d). Some of these genes are known to be required 
for both DNA replication and the DNA checkpoint, and 
mutant strains enter mitosis with incompletely replicated 
DNA [16-20]. The miss E category also contains genes 
required for mitosis, for example, the kinetochore protein 
Nnfl and pericentrin Fcpl. In total, 57 genes were annotated 
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Figure 2. Distribution of cell shape genes, [a] Distribution of 1395 genes with a 
cell shape deletion phenotype among the 11 cell shape categories, (b) Overlap of 
long gene set identified in this study (513 genes, green cirde) with a set of pre- 
viously published genes with a long deletion phenotype (158 genes, red cirde). 
For further details see the electronic supplementary material 1, table S8. (c) Over- 
lap between (i) 521 fission yeast orthologues of human genes with an RNAi cell 
cycle phenotype (blue drde), (ii) 614 genes with a mitotic cell cyde annotation in 
fission yeast (pink drde) and 276 new cell cyde genes from this study (orange 
drde). For further details see the electronic supplementary material 1, tables SB 
and S11. 



to one or more cell cycle process (table 2 and electronic sup- 
plementary material 1, table S6d). We conclude that the miss 
E phenotypic class w^ill be a good starting point to screen for 
further genes required for both DNA replication and the 
DNA checkpoint as well as for genes required during mitosis. 

Another group of previously know^n cell cycle genes exhi- 
bit a wee phenotype, w^ith viable cells dividing at a shorter 
length and a smaller cell volume than WT cells. A partial 
gene deletion library of viable haploid mutants has already 
been screened, and 18 wee mutants identified [21]. In this 
study, we identified by visual examination 25 genes w^ith a 
small cell deletion phenot3rpe; this small set consisted of 
11 non-essential genes and 14 essential genes (see the electronic 
supplementary material 1, table S5m). The 11 non-essential 
genes included nine of the 18 previously identified wee genes 
and were mainly those w^ith a stronger wee deletion phenotype 
[21]. The two additional genes encode a predicted 26S protea- 
some non-ATPase regulatory subunit (SPCC18.17c), w^hich 
has been found (F. Navarro and P. Nurse 2012, personal com- 
munication) to be shorter but also wilder and so divided at a 
WT cell volume, and a mitochondrial inheritance GTPase 
(dmll), w^hich is potentially a new wee gene. The 14 essential 
genes were enriched for tRNA metabolism (see table 2 



Table 1. Summary of phenotype categories. Fourteen broad phenotype 
categories were used for analysis, and the number and dispensability of 
genes in the different categories is shown. Each gene is dassified by a 
single phenotype category based on the most penetrant or strongest 
deletion phenotype. Cells were only described as wild-type (WT) when no 
other phenotype was observed. The dassifier for each gene can be found in 
the electronic supplementary material, table SI, column H and gene 
dispensability in table SI, column I. All genes are induded only once and 
cell phenotype terms from the electronic supplementary material 1, table 
S4 that are not induded as separate phenotypes are subsets of one of 
these categories, for example, a T-shaped mutant is induded in the curved 
category as curved is the most penetrant phenotype for this type of cell 
shape mutant. 
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for summary and the electronic supplementary material 1, 
tables S6m and S7m for details), including genes encoding 
subunits of the RNase P and mitochondrial RNase P, 
w^hich have roles in RNA processing [22], tRNA 2'-0-ribose 
methy transferase, tRNA-specific adenosine deaminase (2/2 
subunits) and the tRNA-specific splicing endonuclease (2/4 
subunits). It is possible that these non-viable smaU mutants, 
like other similar small size viable mutants identified in 
budding yeast that affect grow^th [23], may only indirectly 
affect cell cycle progression. 



3.4. New cell cycle genes 

Our genome-w^ide screen has identified 513 genes w^ith a 
long cell deletion phenotype and thus required for the cell 
cycle in fission yeast. Previously, 158 fission yeast genes 
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Table 5. GO process annotations for 276 novel cell cycle genes. These 276 
genes were not previously known to be involved in the cell cycle in fission 
yeast. Two hundred and thirty genes are annotated to other GO processes 
that have accepted links to the cell cycle in fission yeast and 29 genes had 
a GO process annotation not related to the cell cycle. Only 17 genes were 
of completely unknown function. 
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that generate elongated cells when deleted have been 
reported and annotated in PomBase (http://www.pombase. 
org/) [24]. To validate this qualitative visual approach to 
identify new cell cycle genes, we compared these 158 genes 
(see the electronic supplementary material 1, table S8a) with 
the 513 fission yeast cell cycle genes from this screen (see 
the electronic supplementary material 1, table S5g-i) and 
found that 147 of the 158 genes were also identified in our 
screen (see figure 2h and electronic supplementary material 1, 
table S9). The 366 genes not previously reported as 
elongated when deleted included 90 genes with an existing 
cell cycle GO annotation and 276 genes with no previously 
known cell cycle role (see the electronic supplementary 
material 1, tables S^h and S8c). The majority of these 276 
genes (230/276) are annotated to GO processes which have 
previously been linked to the cell cycle, suggesting that the 
276 genes are true positives and identify new cell cycle 
genes. The 230 new genes involved in cell-cycle-related pro- 
cesses included genes required for ribosome biogenesis, 
splicing and nucleotide metabolism (table 5). For example, 
seven genes {dfrl, adkl, hptl, deal, dutl, dcdl, tmpl) are 
concerned with various nucleotide metabolism pathways. 
Two previously identified cell cycle genes, budding yeast 
CDC8 (tmpl orthologue) [25], and the fission yeast 
cdc22 (ribonucleotide reductase) [26] are also required for 
nucleotide metabolism. Given the cell cycle role of these 
two genes, the other genes identified here may also be impor- 
tant for maintaining the nucleotide levels needed for cell 
cycle progression. 

The remaining 46 genes include 17 genes unstudied in 
any organism (nine of which are conserved in humans). 
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and 29 genes that have existing GO annotations to processes 
or pathways not previously linked to the cell cycle. Of these, 
13 genes are involved in a number of different metabolic 
pathways, including amino acid, carbohydrate and phospho- 
lipid metabolism. We investigated whether any of these 
genes had genetic or physical interactions with genes impli- 
cated in the cell cycle using the BioGRID Interaction 
database [27]. Several genes showed such interactions (see 
the electronic supplementary material 1, table SIO). For 
example, a predicted pyruvate decarboxylase SPAClF8.07c 
interacts with a wee gene zfsl [21,28]. It is possible that 
these 13 metabolic cell cycle genes may act as regulatory 
links between small molecule biosynthesis pathways and 
the cell cycle. 

3.5. Comparison with a human cell cycle gene set 

To examine the overlap between cell cycle genes in fission 
yeast and human, we identified a set of 521 human genes pro- 
posed to be involved in the cell cycle [12], and which have a 
fission yeast orthologue (see the electronic supplementary 
material 1, table SI la and §5.5.3). The 521 fission yeast ortho- 
logues of these human genes were compared with 614 fission 
yeast genes with an existing mitotic cell cycle annotation, 
including all genes so far annotated to the mitotic cell cycle 
either by inference or experiment (see the electronic sup- 
plementary material 1, table Sill? and §5.5.1; http://www. 
pombase.org/). There were 113 genes common to both gene 
sets (see figure 2c and electronic supplementary material 1, 
table Sllc). We also compared the 521 gene set with the 276 
new cell cycle genes from this study and identified a further 
43 genes in common (see figure 2c and electronic supplemen- 
tary material 1, table Slid). Therefore, in total, 156 of the 521 
conserved genes were involved in the cell cycle in both organ- 
isms (29.9%), a similar level to other inter-species comparisons, 
human /worm at 36 per cent and human /fly at 38 per cent 
[12,14,29,30]. Possible reasons why these inter-species com- 
parisons in a variety of studies show such a low overlap are 
considered in §4. 

3.6. Cell cycle checkpoint genes 

To identify new DNA checkpoint genes, we screened 2983 
viable gene deletion mutants for those that failed to block the 
cell cycle in the presence of the ribonucleotide reductase inhibi- 
tor HU [31-33]. We further screened these HU-sensitive 
mutants for those with a cut phenotype [34], where cells fail 
to block cell cycle progress and enter mitosis generating 
chromosome segregation defects. We identified 132 mutants 
sensitive to at least 5 mM HU (see §5.2 and electronic sup- 
plementary material 1, table S12). In the presence of HU, 
deletion mutants of eight genes had greater than 60 per cent 
cells with a cut phenotype (see the electronic supplementary 
material 1, table S12, column O). Of these, deletion mutants 
of husl, mdl, md3, md26, radl7 and rad9, which have been pre- 
viously been shown to be required for establishment of the 
DNA checkpoint, did not elongate prior to entering mitosis 
[33]. Deletion mutants of the remaining two genes {ddhl, 
lem2) initially became elongated but eventually entered mitosis 
and displayed a cut phenot3rpe (figure 3a). Ddbl is necessary 
for stabilizing DNA replication forks and is involved in regulat- 
ing the replication checkpoint kinase Cdsl [35], and a lem2 
mutant has previously been shown to be sensitive to HU 




systematic ID 


new 

checkpoint 
gene 


% cut 
after 

10 h in HU 


SPAC18G6.10 


lem2 


88 


SPAC1805.04 


nupl32 


26 


SPAC19E9.01C 


nup40 


27 


SPAC4F8.01 


did4 


27 


SPAC8C9.17C 


spc34 


23 


SPAC9E9.14 


vps24 


28 


SPBC19C7.02 


ubrl 


24 


SPBC6BL04 


mde4 


28 


SPBP8B7.10C 




22 


SPCC1393.05 


ersl 


20 



Figure 3. New HU-sensitive cut genes, [a) Iem2 deletion mutant cells stained 
with DAPI after growing for 8-10 h in the presence of HU. Examples of anu- 
cleated cells can be seen (white arrow) and cells with unequally segregated 
chromatin (red arrow), (b) New checkpoint genes identified in this study. 



[36]. Mutants of a further 24 genes showed a lower level of 
cut cells (between 20 and 60%) after 10 h in HU (see the elec- 
tronic supplementary material 1, table SI 2, column P). Of 
these, nine genes {nupl32, nup40, did4, spc34, vps24, uhrl, 
mde4, utpl6 and ersl) are newly identified as HU-sensitive 
genes with a cut phenotype. 

In this study, we have identified one new S phase check- 
point gene, lem2, which has a strong cut deletion phenotype, 
and nine genes with a lower penetrance cut deletion pheno- 
type, which may influence maintenance of the S phase 
checkpoint (figure 3b). 

3.7. Cell shape mutants 

Cell shape mutants other than those with a long cell pheno- 
type have been used to identify genes important for 
generating the normal rod-shape of a fission yeast cell. Pre- 
vious studies in fission yeast have identified orb mutants 
that are spherical because cells fail to grow in a polarized 
manner, ban (banana) mutants that have a curved or bent 
cell phenotype because cells no longer orientate the growth 
zones at 180° along the long axis of the cell, and tea (tip 
elongation aberrant) mutants, which form a new growth 
zone in the wrong place often at 90° to the long axis of the 
cell [2,3]. Other cell shape mutants included bottle- or 



skittle-shaped cells and more generally misshapen cells [37]. 
We screened the deletion mutants for shape defects using 
the following seven categories to describe the mutant pheno- 
types (see figure 1 and electronic supplementary material 1, 
table S5d-f,i-l,n for gene lists), (i) rounded, which includes 
the typical orb mutants and also mutants that are more 
rounded than WT but not completely spherical, (ii) stubby, 
which look shorter and wider than WT but are mainly rod- 
shaped. These two categories showed a degree of overlap, 
with some mutants showing both phenotypes. (iii) curved, 
which includes the ban mutants and the tea mutants. 
During vegetative growth, tea mutants have only a low 
level of T-shaped cells with most cells having a curved phe- 
notype. (iv) skittle: one end of cell the is wider than the 
other end. A total of 333 genes (6.9% of total genes) 
showed these specific alterations in cell shape when deleted 
and are thus important for the generation of normal cell 
shape. A further less well-defined group called misshapen 
showed an ill-defined potato-like shape or a mixture of 
other shapes. These fell into three further subgroups: 
(v) viable misshapen mutants (miss V), (vi) viable misshapen 
mutants, which have a weak phenotype {miss weak V) and 
(vii) essential misshapen mutants {miss E). There were 524 
genes with these more general misshapen deletion pheno- 
types. In total, 857/4843 genes (17.7%) showed altered cell 
shape when deleted (table 1) and 668 of these are conserved 
in humans (77.9%). No additional cell shape phenotypes 
were identified compared with earlier work, suggesting that 
there may only be a restricted number of shapes that a fission 
yeast cell can adopt. 

The rounded, stubby and curved sets are all enriched for 
genes implicated in cell polarity and for localization at the 
cell tip (tables 2 and 3), rounded and stubby sets for cell wall 
organization, and the stubby set for cytokinesis and actin cyto- 
skeleton organization (table 2, footnote k). All 14 genes in the 
curved set annotated to the cytoskeleton organization category 
are involved in microtubule cytoskeleton organization (see 
table 2, footnote 1) and are also enriched at the cell tip (table 3). 
Therefore, we predict that unknown genes with a curved 
deletion phenotype when deleted are likely to be involved in 
microtubule-related processes. Similarly, the stubby pheno- 
type is likely to be associated with genes that affect actin 
processes. Genes which generate a skittle phenotype when 
deleted were enriched for mitochondrial organization; 
nearly, 50 per cent (118/241) of the total genes annotated to 
mitochondrion organization were in the skittle category. We 
found that 19 genes were required for mitochondrial tRNA 
metabolism (table 2) and 61 genes for the mitochondrial ribo- 
some, suggesting that mitochondrial translation underlies the 
skittle phenotype. The miss E category was enriched for genes 
required for lipid metabolism (35 genes), 15 of which are 
involved in glycosylphosphatidylinositol (GPI) anchor biosyn- 
thesis (table 2, footnote a; electronic supplementary material 1, 
table S6d). GPI anchor proteins affect cell wall integrity, and 
loss of these proteins can result in a misshapen cell phenotype 
[38]. In higher eukaryotes, GPI anchor proteins have been 
implicated in the sorting of membrane proteins important 
for cell polarization [39], and so, in fission yeast, GPI anchor 
proteins could also be implicated in cell polarization. 

To investigate whether gene deletions that cause cell shape 
changes also have defects in the cytoskeleton, bipolar growth 
pattern or the cell wall, we analysed 54 previously uncharacter- 
ized viable cell shape mutants. These were a subset of 352 



non-essential genes in the miss V, miss weak V, rounded, stubby, 
curved and skittle categories (table 1). We found that a total of 
35 strains had defects in the cytoskeleton and /or the cell wall, 
bipolar growth or cell separation defects and included 26 
mutants with actin or microtubule defects (see figure 4 and 
electronic supplementary material 1, table S13 and electronic 
supplementary material 2 for details). These results indicate 
that further screening of viable deletion mutants, even those 
with a weak phenot3^e, will be a useful approach to identify 
and characterize genes involved in the cytoskeleton. 



4. Discussion 

We have visually screened 4843 gene deletion mutants in fis- 
sion yeast, 95.7 per cent of all protein coding genes, and have 
identified near genome-wide sets of the genes required for 
the cell cycle and cell shape, the first systematic description 
for a eukaryote. The long cell phenotype in fission yeast 
defines cells blocked in cell cycle progression during inter- 
phase or cytokinesis and so is an effective way to identify 
genes required for these stages of the cell cycle. GO enrich- 
ment analysis has shown that the long HP and long Br 
categories were enriched for genes previously identified as 
being required during interphase or cytokinesis respectively. 
The long LP set included genes previously known to be 
required during mitosis and we suggest that these genes 
may represent a subgroup of mitotic genes that have 
additional roles during interphase. Several genes required 
for the cell cycle were found in the miss E gene set suggesting 
that further analysis of this set may also identify new genes 
required during interphase or mitosis. 

We identified 513 cell cycle genes in total, 276 of which 
were not previously known to have a role in the cell cycle. 
Of these new genes, 230 were annotated to GO processes pre- 
viously implicated in the cell cycle, thus identifying new links 
between the cell cycle and these processes. These included 
genes required for nucleocytoplasmic transport, mRNA meta- 
bolic process (specifically splicing), ribosome biogenesis and 
nucleotide metabolism (table 5). There were 46 new cell cycle 
genes not annotated to a process previously associated with 
the cell cycle, 13 of which are involved in small molecule 
metabolism. Frequently, only one or two genes were identified 
for a specific metabolic pathway. For example ectl is predicted 
to encode ethanolamine-phosphate cytidylyltransferase, which 
is rate-limiting for S3mthesis of CDP-ethanolamine, an impor- 
tant step in phospholipid bios3mthesis [40]. We speculate that 
these genes may encode proteins linking different aspects of 
metabolism, such as metabolite levels or flux, to the cell cycle. 

Studies using RNAi in metazoan organisms have ident- 
ified sets of genes required for cell cycle progression. 
However, these overlap only to a limited extent between var- 
ious intra- and inter-species comparisons in a range from 10 
to 38 per cent [12]. Our analysis comparing cell cycle genes 
in fission yeast and human identified 156/521 orthologous 
gene pairs (29.9%) involved in the cell cycle in both organ- 
isms, a similar percentage overlap to that found in other 
inter-species studies. A possible reason for the rather limited 
overlap of cell cycle genes in a wide range of studies may be 
because gene knockdowns using RNAi can result in varying 
levels of gene product and thus more variation in cellular 
phenotype. Our analysis is based on gene deletions that gen- 
erally eliminate the entire gene function, thus reducing 
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Figure 4. Cytological analysis of novel cell shape mutants. Examples of viable mutants from four cell shape phenotype categories analysed for defects in the 
cytoskeleton, growth pattern or cell wall, (a) Wild-type, [b] meu29A.. (c) y^7f9A. (d) spclls.. [e] tlgllS.. DIC, differential interference contrast; CF, calcofluor 
used to stain the cell wall and septum. For details, see the electronic supplementary material 2 and electronic supplementary material 1, table S13. 



phenotypic variability. Inter-species comparisons may also be 
limited because different phenotypes may be used in differ- 
ent organisms to infer a specific cell cycle defect, and these 
are not always directly comparable. A cell cycle defect 
during interphase in fission yeast produces an easily identifi- 
able highly consistent long cell phenotype, which can be used 
to reliably identify interphase cell cycle genes. Many of the 
cell cycle genes identified in fission yeast have a conserved 
cell cycle role in other eukaryotes and, so the long cell cycle 
gene set identified in this study is likely to be useful to 
uncover additional cell cycle genes conserved across species. 

To catalogue genes required to generate and maintain the 
correct cell shape, we identified categories of deletion mutants 
with specific cell shape defects. GO analysis showed that genes 
within these groups could be used to identify genes implicated 
in the actin cytoskeleton (stuhhy), the microtubule cytoskeleton 



(curved) and mitochondrial function (skittle). The skittle cat- 
egory suggests that there is an uncharacterized mechanism 
influencing cell shape, which is affected by defects in mito- 
chondrial organization and translation. The distribution of 
mitochondria in a cell is dependent on microtubules [41], 
and binding of mitochondria to microtubules can moderate 
microtubule dynamics [42], raising the possibility that defec- 
tive mitochondria may affect microtubules thus leading to 
cell shape changes. In humans, a number of diseases includ- 
ing deafness and muscle pathologies are linked to defects in 
mitochondrial protein synthesis [43]. The link we have ident- 
ified in this study between mitochondria and cell shape 
suggests that cell shape changes could be an underlying 
cause of some of these pathologies in humans. We also 
showed that some genes from the miss weak V set, although 
only exhibiting a mild shape change when deleted, have 



cytoskeletal defects. This suggests that genes with this deletion 
phenotype will be a good source of new genes affecting the 
actin and microtubule cytoskeleton. 

A limited range of defined cell shapes were identified in 
the genome-wide gene deletion screen. It appears that only 
a restricted range of cell shapes is possible for the fission 
yeast cell, perhaps reflecting topological constraints in cellu- 
lar organization, related to the cell wall or the cytoskeleton, 
for example. Future comparisons with other organisms that 
have been screened for cell shape defects [10,11,13] will 
help identify the genes and processes required to generate 
or maintain eukaryotic cell shape. 

Our work provides the first near genome-wide sets of 
gene deletions that influence the eukaryotic cell cycle and 
cell shape. This qualitative classification of genes according 
to cell shape phenotypes, based on deletion mutants, pro- 
vides a resource that will be a good starting point for 
further studies in fission yeast and for the identification of 
equivalent gene functions in other eukaryotic organisms. 

5. Material and methods 

5.1. Phenotype analysis of the genome-wide set of 
gene deletions 

The 4843 deletion strain collection used for this analysis 
consists of 4825 strains described by Kim et al. [7] plus 18 
additional gene deletion mutants as shown in the electronic 
supplementary material 1, table S2. Changes in gene dispen- 
sability from Kim et al. [7] for nine reconstructed strains and 
12 re-analysed strains are shown in the electronic supplemen- 
tary materia 1, table S3. All growth conditions and media 
were used as described by Moreno et al. [44], unless other- 
wise stated. Spores were generated as described for gene 
dispensability analysis [7]. All strains were coded, and a 
blind analysis was conducted. Between two and four isolates 
for each heterozygous diploid deletion mutant were indepen- 
dently sporulated. A visual examination of the phenotypes of 
both deletion G418-resistant spores and WT G418-sensitive 
spores following free spore analysis was carried out after 1 
and 2 days following plating on non-selective YES plates at 
25°C and 32°C. The presence of both G418-sensitive and 
-resistant spores allowed a comparison of the deletion 
mutant with WT. Any phenotypic differences from WT that 
could be detected by eye were described as the putative del- 
etion phenotype. After 2 days, colonies were replica plated 
onto YES plates containing G418 (Sigma) at 25 °C and 32 °C 
to confirm the gene deletion phenotype by linkage to the 
G418-resistant phenotype. 

The final deletion phenotype categories for a genome- 
wide set of genes were generated as follows, a GO analysis 
of the long group was compared with a GO analysis of the 
subdivisions long HP, long LP and long Br. These subdivisions 
formed biologically significant subgroups of the long group 
and so these three categories were used for further analysis. 
The same type of GO analysis was used for the rounded 
and stubhy groups as these formed biologically distinct cat- 
egories, although there is also overlap between the 
phenotypes of these two categories. The misshapen group 
was divided into essential and non-essential genes, because 
viable mutants may be more useful for identifying genes 
required for cell shape where as the misshapen phenotype 



observed in the essential group may be less specific for a 
cell shape defect given that the cells are dying or dead. The 
remaining categories WT, spores, germination, skittle, curved 
and small could not be usefully further subdivided by their 
phenotype. 

To estimate the minimal cell length increase detectable, we 
measured the cell length of 34 viable gene deletion mutants 
described as long after visual screening on plates and which 
had not previously been implicated in the cell cycle (see the 
electronic supplementary material 1, table SI column J). We 
could detect cells at least 10 per cent or longer compared 
with WT (approx. 15.6 ixm) and, so a 10 per cent or greater 
increase in cell length was used as the criteria for a long phe- 
notype. The cut-off between high penetrance and low 
penetrance was 30 per cent long cells. This was estimated 
using inviable mutants that formed microcolonies showing a 
mixture of long and WT/short cells. For these mutants, 
30 per cent or less of the cells had a long phenotype. 

To validate this approach, we compared the 513 genes 
from the three long categories with 158 cell cycle genes 
reported in PomBase as long. We identified 147/158 (93%) 
of these genes, suggesting that the remainder of the genes 
in our long category are also likely to be involved in the cell 
cycle. Furthermore, only 46/513 were not annotated to a 
GO category previously linked to the cell cycle. 

CeUs showing different cell shape defects were photo- 
graphed using a Zeiss Axioskop microscope with a CF plan 
X50/0.55 objective and a Panasonic DMC-LX2 camera. 
Spores from representative strains were plated on to YES 
solid medium and allowed to germinate or form small colonies 
before being photographed. 

5.2. Screen for new DNA checkpoint genes 

The growth of 2983 viable deletion mutants from Bioneer ver- 
sion 1 were screened on YES agar plates for 24-48 h either in 
the presence of 2.75 or 5.5 mM HU or without HU and scored 
on a scale of strong (+++), medium (++), weak (+) or no sen- 
sitivity, depending on their ability to grow on different HU 
concentrations compared with no HU (see the electronic sup- 
plementary material 1, table S12). To check whether the 132 
HU-sensitive mutants were also involved in the DNA check- 
point preventing mitosis, cells were grown in liquid cultures 
with 11 mM HU and screened for a cut phenotype using 
4',6-diamidino-2-phenylindole (DAPI) to visualize the nucleus. 

5.3. Cell length measurements 

CeUs were grown to mid-exponential growth (2 x 10^ to 1 x 
lO^ceUs ml"^) in YES liquid medium at 32°C (or 25°C where 
appropriate) and photographed using a Zeiss Axioplan micro- 
scope with 100 X objective and a COHU CCD camera. Cell 
lengths of 30 septated cells were measured using ImageJ. 

5.4. Phenotypic analysis and cytoskeleton analysis of 
viable shape mutants 

For the initial characterization, cells were grown at 18°C, 25.5°C, 
29°C and 34°C on minimal and YES medium plates, and ceU 
morphology analysed by differential interference contrast 
(DIC) microscopy. For further characterization, strains showing 
morphological defects were grown in liquid rich medium at 



25°C to mid-log phase or in the conditions at which each strain 
showed the strongest phenot3rpe by DIC microscopy. Septa 
and pattern of cell growth was visualized with 35 |jig ml~^ 
calcofluor staining (fluorescent brightener; Sigma). For WT 
cells, the septation index was 15 per cent (n = 500 cells) 
and 30.1 per cent of cells showed monopolar growth (n = 300 
cells). Nuclei were visualized with 0.2|jLgml~^ DAPI 
(Sigma) or 100 |jim ml~^ IP (Sigma) staining. Actin staining 
was as described by Pelham & Chang [45] using AlexaFluor 
488-phaUoidin (Molecular Probes). For anti-tubulin immuno- 
fluorescence, cells were fixed in methanol at — 80°C and 
further processed as described by Hagan & Hyams [46]. Primary 
antibodies were anti-tubuHn ((TAT-1; 1 : 80 dilution) followed by 
Alexa 488 goat anti-mouse secondary antibody (Molecular 
Probes). Microscopy was performed at 23-25°C, either with an 
Axioplan 2 microscope (Carl Zeiss, Inc.) equipped with a Cool- 
snapHQ camera (Roper Scientific) or with a Leica TCS SL 
confocal microscope. Data were acquired using the 100 x 
objective taking seven z-sections with 0.5 ixm spacing. 

5.5. Bioinformatics analysis 

5.5.1. Identification of genes already implicated in mitotic cell 
cycle processes 

To identify genes involved in the mitotic cell cycle in fission 
yeast, we used fission yeast GO data from 26 September 
2011 (http://www.pombase.org/), and selected the set of 
protein coding genes annotated to: 

00:0000278 mitotic cell cycle, 

00:0000910 cytokinesis, 

00:0006261 DNA-dependent DNA replication, and 
00:0000075 cell cycle checkpoint. 

Minor adjustments were made to this dataset to remove 
three known false positives (SPAC343.17c, SPBC19F8.02, 
SPAC5D6.08c) and add three known false negatives 
(SPAC23H4.11C, SPAC26A3.03C SPAC23H4.18c). The 
complete list is provided in the electronic supplementary 
material 1, table Sill?. 



5.5.2. Gene Ontology enrichment analysis 

OO enrichments were performed using OO term finder 
(http:// go.princeton.edu/ cgi-bin/ OOTermFinder with ontol- 
ogy and annotations from 26 September 2011). Threshold 
pvalues of 0.001, 0.01 and 0.1 were used to identify specific 
enrichments. Bonferroni correction was used. 

OO slim categories presented in tables 2, 3 and 4 refer to 
the OO IDs in the electronic supplementary material 1, table 
S14. All OO terms and p-values for each phenotype set are 
provided in the electronic supplementary material 1, tables 
S6 and S7. Oenes where number of annotations = 1 were 
obtained using OO term Mapper (http://go.princeton.edu/ 
cgi-bin/ OOTermMapper). 

The background set was the 4843 gene set used this study. 



5.5.3. Comparison with human cell cycle genes 

A list of 1351 human cell cycle genes was extracted from a 
study by Kittler et al. [12] and mapped to current Ensembl 
IDs (79 of the human identifiers had been retired and were 
no longer linked to extant genes). The remainder were 
mapped to 521 fission yeast orthologues using Ensembl 
Compara 10/11/2011 (http://genome.cshlp.org/content/ 
19/2/327.1ong). 
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